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ABSTRACT 

A 

Broadband  seismic  waveform  data  recorded  in  the  distance  range 
of  1  to  11  km  from  two  underground  explosions  detonated  at  Pahute 
Mesa  of  the  Nevada  Test  Site  have  been  used  to  estimate  the  source 
time  functions  for  the  second-order  moment  tensors  of  these  events. 
Data  were  available  from  8  three-component  stations  for  the  explosion 
Harzer  and  from  11  stations  for  the  explosion  Chancellor.  Stable 
moment  tensor  estimates  with  a  good  signal-to-noise  ratio  were  obtained 
in  the  frequency  range  of  0.2  to  5.0  Hz.  The  moment  tensors  are  dom¬ 
inated  by  the  diagonal  elements,  but  a  significant  asymmetry  is  present, 
with  vertical  force  couples  having  a  different  time  dependence  than  the 
horizontal  force  couples.  When  the  isotropic  parts  of  the  moment  ten¬ 
sors  are  interpreted  as  reduced  displacement  potentials  the  results  are  in 
reasonable  agreement  with  theoretical  and  empirical  models,  although 
the  parameters  which  are  estimated  with  these  models  show  systematic 
differences  from  published  scaling  relationships.  These  differences  may 
be  related  to  the  fact  that  the  explosions  were  detonated  in  unsaturated 
materials  or  to  the  effects  of  spall  in  surficial  layers  above  the  source. 
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1.  INTRODUCTION 


Although  it  would  appear  to  be  a  relatively  simple  problem,  there  is  still  much 
that  is  not  understood  about  the  generation  of  elastic  waves  by  large  underground 
explosions.  In  general  the  seismograms  produced  by  such  explosions  are  considerably 
more  complicated  than  what  one  would  expect  for  a  simple  model,  such  as  a  pressure 
pulse  applied  to  the  interior  of  a  spherical  cavity  in  a  homogeneous  halfspace.  Pan  of 
the  difficulty  may  be  related  to  the  fact  that  the  explosions  are  detonated  in  media 
which  are  heterogeneous  in  terms  of  both  physical  properties  and  preexisting  stresses 
and  there  is  dynamic  interaction  between  the  explosion  and  these  heterogeneities, 
including  the  free  surface.  Also  contributing  to  the  difficulty  is  the  media  which  lies 
between  the  source  region  and  the  receiver  where  the  waves  are  recorded,  for  this 
media  is  generally  heterogeneous  and  the  associated  wave  propagation  effects  can  add 
considerable  complexity  to  the  seismograms.  It  would  appear  that  the  first  step  in 
striving  to  obtain  a  better  understanding  of  the  explosions  source  process  is  to  remove 
as  much  as  possible  of  the  wave  propagation  effects  due  to  the  intervening  media 
between  source  and  receiver. 

An  objective  of  the  present  study  is  to  obtain  unobstructed  estimates  of  the  source 
processes  of  large  underground  explosions.  Such  estimates  are  obtained  through  the 
method  of  moment  tensor  inversion  in  which  calculated  Green  functions  are  used  to 
remove  propagation  effects  and  the  source  is  characterized  in  terms  of  the  second-order 
moment  tensor,  a  fairly  general  mathematical  representation  of  spatially  localized 
seismic  sources.  Having  obtained  estimates  of  the  second-order  moment  tensor,  these 
can  be  interpreted  in  terms  of  the  physical  processes  in  the  source  region,  which  is  a 
second  objective  of  the  study. 

The  existence  of  numerous  waveform  data,  previous  studies  in  the  same  area,  and 
an  accurate  velocity  model  all  combine  to  make  Pahute  Mesa  a  convenient  site  for  the 
present  study.  Pahute  Mesa,  a  portion  of  the  Nevada  Test  Site,  has  been  the  location 
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of  numerous  underground  nuclear  explosions  in  the  past  and  continues  to  be  an  active 
location  for  the  detonation  of  buried  explosions.  A  large  amount  of  ground  motion 
data  have  been  collected  from  these  explosions,  with  numerous  recordings  at  local, 
regional,  and  teleseismic  distances.  These  data  have  been  the  subject  of  several 
scientific  investigations  which  provide  useful  background  information  and  also  provide 
results  which  can  be  compared  with  those  of  the  present  study.  Some  of  these  previ¬ 
ous  studies  of  particular  interest  include  Helmberger  and  Hadley  (1981),  Stump  and 
Johnson  (1984),  Barker  et  al.  (1985a,  1985b),  and  Patton  (1985).  The  results  of  the 
moment  tensor  inversion  method  depend  critically  upon  the  availability  of  accurate 
Green  functions  and  this  requires  that  an  accurate  velocity  model  be  available.  A 
recent  study  by  Leonard  and  Johnson  (1987)  has  provided  such  a  model  for  Pahute 
Mesa. 

The  data  used  in  the  present  study  were  all  recorded  within  11  km  of  the  explo¬ 
sion  epicenters.  There  are  certain  advantages  of  using  data  from  these  close  distances. 
First,  many  propagation  effects,  both  elastic  and  inelastic,  can  be  minimized  by 
remaining  as  close  as  possible  to  the  source.  This  can  be  important  when  trying  to 
recover  source  information  at  the  high  frequencies  where  attenuation  and  scattering  can 
have  significant  effects.  Second,  at  these  distances  the  seismograms  contain  P  and  S 
phases  which  left  the  source  travelling  in  both  upward  and  downward  directions  as 
well  as  both  Rayleigh  and  Love  surface  waves.  Thus  the  focal  sphere  is  sampled 
fairly  well,  with  the  major  gap  in  coverage  being  the  P  waves  that  leave  the  source 
region  travelling  downward  at  small  angles  of  incidence  and  emerge  at  teleseismic  dis¬ 
tances. 
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2.  DATA 


The  explosions  Harzer  and  Chancellor  were  both  detonated  in  Area  19  of  Pahute 
Mesa  at  the  Nevada  Test  Site.  The  parameters  of  these  events  are  listed  in  Table  1. 
Note  that  the  two  events  are  only  separated  by  4.3  km  and  the  shot  environments  are 
similar. 

Ground  motion  recordings  at  near  distances  to  these  events  were  obtained  with 
triggered  digital  event  recorders.  Each  instrument  package  consisted  of  three  force- 
balance  accelerometers  oriented  to  record  radial,  transverse,  and  vertical  ground 
motion.  The  data  were  recorded  at  200  samples  per  second  with  a  resolution  of  12 
bits  per  sample  and  the  complete  system  was  flat  to  acceleration  between  frequencies 
of  0.2  and  50  Hz.  The  signal-to-noise  ratio  generally  attains  a  maximum  of  about 
1000  in  the  frequency  band  of  1  to  10  Hz,  but  decreases  at  both  lower  and  higher  fre¬ 
quencies.  Because  the  recording  instruments  were  accelerometers,  this  ratio  falls  to  a 
value  of  1  at  a  frequency  of  about  0.2  Hz,  which  means  that  periods  longer  than  5  sec 
are  not  well  resolved  in  these  data. 

The  locations  of  the  two  events  and  recording  sites  are  shown  in  Figure  1.  Data 
were  obtained  at  8  different  stations  for  Harzer  and  at  1 1  stations  for  Chancellor.  The 
station  distances,  azimuths,  and  elevations  are  listed  in  Tables  2  and  3  along  with  the 
maximum  accelerations.  Both  experiments  were  conducted  within  the  confines  of 
Silent  Canyon  Caldera,  a  thick  section  of  Cenozoic  volcanic  sediments. 

The  process  of  estimating  the  moment  tensor  from  far-field  acceleration  record¬ 
ings  requires  that  three  numerical  integrations  be  performed.  To  ease  this  process 
somewhat,  the  acceleration  data  were  first  converted  to  velocity  data  by  applying  one 
numerical  integration  before  beginning  the  inversion  procedure.  The  knowledge  that 
the  long-time  velocity  level  must  return  to  zero  prevents  base-line  uncertainties  from 
causing  any  serious  problems  with  this  integration. 


Previous  studies  employing  similar  ground  motion  data  from  Pahute  Mesa 
(McLaughlin  et  al.,  1983)  have  shown  that  spatial  coherence  of  the  waveforms  over 
distances  of  a  few  hundred  meters  decreases  significantly  for  frequencies  above  5  Hz. 
Because  of  this  and  the  fact  that  it  is  difficult  and  expensive  to  calculate  Green  func¬ 
tions  at  high  frequencies,  the  data  were  low-pass  filtered  with  a  two-pole  comer  at  10 
Hz  and  then  decimated  by  a  factor  of  4  so  as  to  reduce  the  Nyquist  frequency  from 
100  Hz  to  25  Hz. 

The  first  ten  seconds  of  velocity  data  for  the  Harzer  and  Chancellor  events  are 
shown  in  Figures  2  and  3,  respectively,  with  the  stations  arranged  in  order  of  increas¬ 
ing  epicentral  distance.  The  records  are  fairly  simple  and  of  short  duration  at  short 
distances  but  both  the  complexity  of  the  records  and  their  duration  increase  with  dis¬ 
tance.  Inspection  of  the  records  at  the  closest  distances  is  sufficient  to  establish  that 
the  duration  of  the  source  time  function  is  less  than  5  sec  for  the  frequency  band  con¬ 
tained  in  these  records.  The  first  motion  on  the  radial  and  vertical  components  are  in 
the  same  direction  at  all  stations,  away  and  up,  which  is  expected  for  a  simple  explo¬ 
sive  source.  More  surprising  are  the  transverse  components,  which  have  maximum 
amplitudes  comparable  to  those  on  the  radial  and  vertical  components  and  at  some  sta¬ 
tions  have  impulsive  first  arrivals  at  the  time  of  the  direct  P  wave,  which  is  not 
expected  for  a  simple  explosive  source.  The  first  motions  on  the  transverse  com¬ 
ponents  have  mixed  polarities  (Tables  2  and  3)  but  a  systematic  lobed  radiation  pattern 
is  not  apparent.  Although  the  scatter  is  considerable,  the  amplitudes  on  all  three  com¬ 
ponents  decrease  roughly  as  the  inverse  epicentral  distance  to  the  first  power. 

There  is  considerable  evidence  present  in  the  seismograms  of  Figures  2  and  3 
which  indicates  that  either  the  source  radiation  pattern  changes  rapidly  with  azimuth  or 
that  local  effects  near  the  recording  sites  can  significantly  influence  the  waveforms. 
For  instance,  stations  HI  and  H3  are  separated  by  only  0.1  km  in  range  and  64  deg  in 
azimuth  yet  have  markedly  different  waveforms  on  all  three  components.  Stations  H6 
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and  H8  show  the  same  type  of  disparity.  A  comparison  of  the  records  from  H5  and 
C4  tends  to  favor  the  hypothesis  of  a  recording-site  effect,  because  these  stations, 
which  were  about  the  same  distances  from  the  two  explosions  and  located  at  sites 
within  1  km  of  each  other,  show  strong  similarities  in  waveforms. 

Station  H3  had  a  partial  malfunction  in  the  field  which  caused  the  signal  recorded 
on  magnetic  tape  to  have  a  very  small  amplitude.  Through  special  steps  taken  in  the 
laboratory  it  was  possible  to  recover  the  records  shown  in  Figure  2,  which  appear  to 
be  acceptable  in  the  time  domain.  However,  when  these  records  are  inspected  in  the 
frequency  domain  they  show  a  signal-to-noise  ratio  at  low  frequencies  which  is  an 
order  of  magnitude  worse  than  the  records  recorded  at  other  stations.  This  could  have 
a  strong  effect  upon  a  moment  tensor  inversion  which  is  performed  in  the  frequency 
domain,  so  the  records  from  station  H3  were  not  included  in  any  of  the  moment  tensor 
inversions. 


3.  RESULTS 


3.1.  Moment  Tensor  Inversions 

The  procedure  used  to  estimate  the  moment  tensors  of  the  two  explosions  follows 
that  outlined  by  Stump  and  Johnson  (1977,  1984).  The  basic  assumption  involved  in 
this  procedure  is  that  the  source  can  be  represented  by  a  second-order  force  moment 
tensor  and  that  higher  order  terms  can  be  ignored.  This  would  appear  to  be  a  valid 
assumption  for  explosions,  since  at  least  the  initial  size  of  the  source  is  small  com¬ 
pared  to  the  wavelengths  of  interest.  The  inversions  are  performed  in  the  frequency 
domain  where  the  relation  between  the  waveform  data  and  the  moment  tensor  elements 
is  linear.  No  assumptions  are  invoked  with  regard  to  any  relationships  between  the 
amplitudes  of  the  various  moment  tensor  elements  or  their  time  functions.  The  inver¬ 
sions  were  achieved  through  a  least-squares  minimization  of  the  difference  between 
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observed  and  predicted  seismograms.  The  equations  of  condition  were  solved  by 
singular  value  decomposition. 

Critical  to  the  moment  tensor  inversions  are  the  Green  functions  which  are  used 
to  remove  propagation  effects  from  the  seismograms.  The  accuracy  of  the  final  results 
is  directly  related  to  the  accuracy  of  the  Green  functions  which  depends  upon  having 
available  a  good  model  of  the  velocity  and  density  structure.  Leonard  and  Johnson 
(1987)  have  recently  produced  estimates  of  the  average  P  and  S-wave  velocity  struc¬ 
ture  in  the  shallow  crust  of  Silent  Canyon  Caldera.  These  results  differ  from  previous 
ones  primarily  in  that  velocities  are  modeled  as  continuous  functions  of  depth  rather 
than  a  series  of  constant  velocity  layers.  The  density  measurements  from  borehole 
sampling  plus  the  results  in  Snyder  and  Carr  (1984)  were  used  to  associate  a  density 
model  with  the  velocity  models.  The  final  results  for  the  upper  4  km  of  the  Silent 
Canyon  Caldera  are  shown  in  Figure  4.  Also  shown  in  Figure  4  are  estimates  of  ane- 
lastic  properties  which  were  used  in  the  calculations.  No  direct  measurements  of  these 
properties  were  available,  so  the  estimates  were  based  primarily  on  previous  experi¬ 
ence  with  shallow  deposits  (Johnson  and  Silva,  1981).  The  quality  factor  for  S  waves 
ranges  from  a  value  of  30  at  the  surface  to  a  value  of  75  at  a  depth  of  4  km,  while  the 
quality  factor  for  P  waves  rangt  „  from  50  to  150  over  the  same  range.  Note  that  these 
quality  factors  should  be  interpreted  to  include  both  the  effects  of  intrinsic  attenuation 
and  scattering,  as  the  latter  effect  appears  to  be  significant  at  high  frequencies  in  the 
Silent  Canyon  Caldera  (McLaughlin  et  al.,  1983). 

It  should  be  noted  that  the  Green  functions  of  this  study  were  calculated  using  a 
single  one-dimensional  model  of  the  shallow  structure  at  Pahute  Mesa.  This  does  not 
imply  that  such  a  model  is  completely  appropriate  for  this  region.  In  fact,  there  exists 
considerable  geological  and  seismological  evidence  to  indicate  that  appreciable  lateral 
heterogeneity  may  be  present  in  the  Silent  Canyon  Caldera  (Leonard  and  Johnson, 
1987).  The  variations  in  waveforms  in  Figures  2  and  3  which  were  discussed  above 


support  this  inference.  Thus  an  argument  can  be  made  for  the  use  of  a  three- 
dimensional  model  for  the  calculation  of  the  Green  functions  or  at  least  for  the  use  of 
a  different  suitably  averaged  one-dimensional  structure  for  each  source-receiver  combi¬ 
nation.  However,  the  information  necessary  to  construct  such  detailed  velocity  and 
density  models  of  the  three-dimensional  structure  are  not  available  and  not  likely  to 
become  available  in  the  near  future.  In  such  a  situation  one  must  appeal  to  some  type 
of  statistical  averaging  approach  in  order  to  mitigate  the  effects  of  the  lateral  hetero¬ 
geneity.  By  including  a  sufficient  number  of  seismograms  recorded  over  a  wide  range 
of  distances  and  azimuths  and  assuming  that  waveform  effects  due  to  the  lateral 
heterogeneity  will  contribute  to  the  moment  tensor  inversion  in  an  incoherent  manner 
while  waveform  effects  due  to  the  average  one-dimensional  structure  will  contribute  in 
a  coherent  manner,  one  reasons  that  inaccuracies  in  the  Green  functions  due  to  lateral 
heterogeneity  may  not  all  be  translated  directly  into  inaccuracies  in  the  estimated 
moment  tensors.  While  such  an  approach  is  definitely  not  as  good  as  having  more 
accurate  Green  functions,  it  does  appear  to  give  reasonable  results.  In  effect,  one  is 
trading  off  inaccuracies  in  the  Green  functions  against  the  number  of  seismograms 
used  in  the  inversion  in  the  hope  that  effects  due  to  the  lateral  heterogeneity  (or  any 
other  inaccuracies  in  the  problem)  will  cancel  themselves  out  in  a  random  manner. 
While  in  principle  the  inverse  problem  of  estimating  the  6  elements  of  the  moment 
tensor  is  completely  determined  given  6  seismograms  (Stump  and  Johnson,  1977),  in 
applications  with  real  data  it  has  been  found  that  about  15  seismograms  are  generally 
required  to  insure  stable  results.  Thus,  in  this  study  21  seismograms  are  used  to  esti¬ 
mate  the  moment  tensor  for  Harzer  and  33  are  used  for  Chancellor. 

The  Green  functions  were  calculated  by  fairly  standard  methods  that  follow  the 
general  outline  of  Kind  (1978).  The  stress-displacement  boundary  conditions  were 
solved  in  the  frequency-wavenumber  domain  using  the  general  method  of  Haskell 
(1953)  with  suitable  high  frequency  modifications  as  described  by  Harkrider  (1964). 
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The  inverse  frequency  and  wavenumber  transforms  were  performed  with  the  trapezoid 
rule  using  a  sampling  rate  in  frequency  governed  by  the  Nyquist  frequency  of  25  Hz 
and  a  sampling  in  wavenumber  governed  by  the  maximum  distance  and  fastest  veloci¬ 
ties  (Alekseev  and  Mikhailenko,  1980).  Attenuation  was  included  according  to  the 
method  of  Silva  (1976).  The  continuous  velocities  and  density  of  the  model  were 
approximated  by  constant  property  layers  of  thickness  50  m  in  the  upper  1.5  km  and 
100  m  at  greater  depths. 

The  moment  tensor  inversions  were  obtained  using  standard  singular  value 
decomposition  techniques  following  the  method  described  by  Stump  and  Johnson 
(1977).  Estimates  of  the  moment  tensor  are  obtained  in  the  frequency  domain  but  are 
also  transformed  back  to  the  time  domain  for  further  analysis  there.  The  inversions 
were  actually  performed  for  the  moment-rate  tensors  and  these  were  then  numerically 
integrated  to  obtain  the  moment  tensors.  In  order  to  provide  sufficiently  dense  sam¬ 
pling  in  frequency  and  also  avoid  possible  truncation  and  windowing  effects,  20  sec  of 
data  were  included  in  the  inversion  procedure.  However,  because  the  transient  signals 
on  the  moment  tensors  were  confined  to  the  first  couple  of  sec  and  because,  as  men¬ 
tioned  above,  the  signal-to-noise  ratio  is  suspect  at  long  periods,  only  the  first  five  sec 
of  the  moment  tensors  have  been  analysed  in  the  time  domain  part  of  this  study. 

The  net  result  of  the  inversion  procedure  is  a  symmetric  second-order  tensor  in 
which  the  6  elements  all  have  independent  time  histories.  The  moment  tensor  M  can 
be  treated  as  a  3  by  3  matrix  with  elements  M n,  Ml2,  —  M33,  where  the  coordinate 
system  is  chosen  so  that  1  denotes  North,  2  denotes  East,  and  3  denotes  down.  To  aid 
in  the  physical  interpretation  of  this  tensor  it  is  helpful  to  decompose  it  into  an  isotro¬ 
pic  part  Mj  and  deviatoric  part  M0,  with  the  definitions 


M  =  Mi  I  +  Md 


where  I  is  a  diagonal  unity  matrix, 
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_  trace  M  _  (M n  +  M 22  +  M33) 
3  3 


Md  =  M  -  M,  I 


The  isotropic  part  can  be  associated  with  the  symmetrical  parts  of  the  source  in  which 
there  is  only  a  change  in  volume,  and  this  is  the  type  of  source  expected  for  a  simple 
explosion.  The  deviatoric  part  can  be  associated  with  asymmetrical  effects  of  the 
source,  such  as  are  expected  for  dislocation  models  of  faulting  or  asymmetrical  parts 
of  an  explosion.  The  interpretation  is  further  aided  by  analyzing  the  deviatoric  part  in 
terms  of  its  eigenvalues  and  eigenvectors.  These  can  be  thought  of  as  the  net  effective 
magnitudes  and  associated  directions,  respectively,  of  the  principal  forces  acting  in  the 
source  region.  Note  that  this  decomposition  in  terms  of  isotropic  parts,  deviatoric 
parts,  and  principal  forces  is  performed  independently  at  each  instant  of  time. 

Before  proceeding  to  a  discussion  of  the  estimates  of  the  moment  tensors  which 
emerged  from  the  inversion  process,  it  is  worth  considering  a  few  details  of  the  pro¬ 
cess  itself.  With  a  total  of  21  seismograms  for  the  Harzer  event  and  33  for  Chancel¬ 
lor,  the  inversions  were  quite  stable  in  a  numerical  sense.  For  instance,  the  ratio  of 
largest  to  smallest  eigenvalue  in  the  singular-value  decomposition,  which  can  be  used 
as  an  estimate  of  the  condition  number,  ranged  between  4  and  27  for  Harzer  and 
between  4  and  16  for  Chancellor.  The  reliability  of  the  inversion  results  depends  on 
frequency  and  this  is  illustrated  in  Figure  5.  Here  we  see  typical  results  for  the  esti¬ 
mates  of  the  moduli  of  the  moment  tensor  elements  and  estimates  of  their  standard 
errors,  which  were  calculated  by  assuming  that  the  noise  in  the  input  data  was  white. 
Note  that  the  ratio  of  the  moment  tensor  estimates  to  their  standard  errors  are  a  max¬ 
imum  in  the  range  of  1  to  5  Hz  and  decreases  at  both  lower  and  higher  frequencies. 
This  is  interpreted  to  mean  that  the  reliability  of  the  moment  tensor  elements  in  best  in 
the  1  to  5  Hz  range  but  deteriorates  at  frequencies  outside  this  range. 
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As  can  be  seen  in  Figure  5,  the  spectral  patterns  apparent  in  the  moment  tensors 
are  quite  similar  for  the  two  explosions,  with  Harzer  having  in  general  slightly  higher 
spectral  amplitudes.  The  M  n  elements  have  reasonably  well  defined  low  frequency 
levels,  comer  frequencies  near  2  Hz,  and  high  frequency  slopes  between  -2  and  -3. 
The  M22  elements,  which  are  not  shown,  are  similar  to  this.  The  M 33  elements  differ 
somewhat  in  that  they  have  an  additional  low  frequency  peak  in  the  0.3-0.4  Hz  range. 
The  average  deviatoric  elements  have  spectral  amplitudes  that  are  about  5  times 
smaller  than  the  isotropic  elements,  have  comer  frequencies  near  0.4  Hz,  and  have 
high-frequency  slopes  of  about  -2.  The  significance  of  these  spectral  differences  in  the 
elements  of  the  moment  tensors  will  be  discussed  more  completely  later  in  the  paper. 

Another  way  of  evaluating  the  moment  tensor  elements  is  to  consider  how  well 
they  explain  the  data.  Such  results  are  shown  in  Figures  6  where  the  observed  seismo¬ 
grams  are  compared  with  predicted  seismograms  calculated  by  convolving  the  Green 
functions  with  the  estimated  moment  tensors  for  two  of  the  Chancellor  stations.  The 
seismograms  from  C3  at  an  epicentral  distance  of  2.1  km  are  quite  compact  in  time 
with  the  direct  P  and  S  phases  having  the  largest  amplitudes,  and  the  estimated 
moment  tensor  for  Chancellor  does  a  good  job  of  explaining  all  of  the  major  features 
on  the  seismogram.  The  correlation  coefficients  between  the  observed  and  calculated 
seismograms  for  the  radial,  transverse,  and  vertical  components  are,  respectively,  0.91, 
0.86,  and  0.81.  The  seismograms  for  station  C7  at  an  epicentral  distance  of  9.8  km  in 
Figure  6  are  typical  of  the  results  at  larger  distances  where  the  seismograms  are  more 
extended  in  time  with  both  surface  waves  and  scattered  body  waves  undoubtedly 
accounting  for  part  of  the  energy.  The  estimated  moment  tensor  does  a  reasonably 
good  job  of  explaining  the  first  4  sec  of  the  radial  and  vertical  components  where  the 
direct  body  waves  are  expected  to  arrive,  but  does  poorly  for  later  times  on  these  com¬ 
ponents  and  does  poorly  for  all  times  on  the  transverse  components.  The  correlation 
coefficients  between  the  observed  and  calculated  seismograms  are  0.41,  0.26,  and  0.24 


for,  respectively,  the  radial,  transverse,  and  vertical  components.  Considering  all  of 
the  Harzer  seismograms,  the  average  correlation  coefficient  is  0.63  ±  0.22  for  the 
radial  components,  0.40  ±  0.30  for  the  transverse  components,  and  0.44  ±  0.22  for  the 
vertical  components.  Similar  average  correlation  coefficients  for  the  Chancellor 
seismograms  are  0.48  ±  0.30,  0.23  ±  0.40,  and  0.37  ±  0.32.  While  correlation 
coefficients  are  insensitive  to  absolute  amplitudes.  Figure  6  shows  that  the  observed 
and  calculated  amplitudes  are  in  reasonable  agreement  on  the  radial  and  vertical  com¬ 
ponents,  with  a  slight  tendency  for  the  calculations  to  underestimate  at  short  distances 
and  overestimate  at  larger  distances.  The  transverse  component  is  consistently 
underestimated  at  all  distances.  Thus,  on  the  average,  the  radial  components  are  fit 
best,  the  transverse  components  are  fit  worst,  and  Harzer  is  fit  better  than  Chancellor. 

Based  on  these  results  showing  the  ability  of  the  moment  tensors  to  fit  the 
observed  data,  a  plausible  interpretation  is  that  the  inversion  process  is  successful  in 
fitting  the  early  parts  of  the  seismograms  where  the  arrivals  are  direct  body  waves 
which  are  likely  to  be  consistent  with  the  Green  functions,  and  the  process  is  also  suc¬ 
cessful  in  ignoring  some  of  the  large  amplitudes  during  the  later  parts  of  the  seismo¬ 
gram  which  are  not  consistent  with  the  Green  functions  and  are  most  likely  due  to 
scattering  and  muti-pathing.  Experience  indicates  that  this  type  of  success  is  depen¬ 
dent  on  having  a  large  number  of  seismograms  so  that  the  problem  is  highly  over¬ 
determined.  As  the  number  of  seismograms  used  in  the  inversion  process  decreases, 
the  likelihood  increases  that  a  random  phase  on  a  particular  seismogram  will  be  inter¬ 
preted  as  part  of  the  source  time  function.  It  follows  that  the  ability  to  fit  the  observed 
seismograms  may  actually  improve  as  the  number  of  seismogams  used  in  the  inversion 
process  is  decreased.  An  implication  of  this  interpretation  is  that  good  agreement 
between  the  observed  and  fitted  seismograms  is  not  necessarily  an  indication  that  the 
moment  tensor  estimates  are  accurate. 
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While  the  moment  tensors  are  actually  estimated  in  the  frequency  domain,  both 
amplitude  and  phase  arc  available  so  that  the  results  can  be  easily  transformed  to  the 
time  domain,  where  the  interpretation  is  often  more  straightforward.  The  first  five 
seconds  of  the  moment  tensors  which  were  estimated  for  the  explosions  Harzer  and 
Chancellor  are  shown  in  Figures  7  and  8,  respectively.  Several  general  observations 
apply  to  the  results  for  both  Harzer  and  Chancellor.  The  diagonal  elements  of  the 
moment  tensors  dominate  and  are  qualitatively  similar  in  both  their  magnitudes  and 
time  histories.  They  are  basically  slightly  acausal  step  functions  with  rise  times  of  less 
than  1  sec.  There  is  discemable  overshoot  on  all  of  the  diagonal  elements  but  this 
effect  is  most  developed  on  the  M33  elements.  On  the  other  hand,  the  off-diagonal 
elements  are  smaller  in  magnitude  than  the  diagonal  elements  and  more  variable  in 
their  time  histories.  In  general,  they  are  more  emergent  and  appear  to  have  less  high- 
frequency  content  than  the  diagonal  elements.  Their  onset  is  generally  delayed  with 
respect  to  the  onset  of  the  diagonal  elements. 

The  differences  between  the  isotropic  and  deviatoric  parts  of  the  moment  tensors 
displayed  in  Figures  7  and  8  are  consistent  with  the  spectral  data  of  Figure  5.  There  is 
a  suggestion  in  these  spectral  data  that  two  separate  functions  might  be  present  in  the 
sources.  The  first  is  a  low  frequency  function  with  a  corner  frequency  near  0.4  Hz  and 
a  spectral  slope  less  steep  than  -2  which  is  present  in  both  the  isotropic  and  deviatoric 
parts.  The  second  is  a  higher  frequency  function  with  a  comer  frequency  near  2.0  Hz 
and  a  spectral  slope  more  steep  than  -2  which  is  present  only  in  the  isotropic  part  of 
the  moment  tensors.  Further  inspection  of  the  spectra  of  the  individual  elements  of  the 
moment  tensor  indicates  that  the  relative  amplitude  of  the  low  frequency  function  is 
larger  on  the  M33  elements  than  on  the  Mu  and  M2 2  elements.  These  results  are 
entirely  consistent  with  the  time -domain  plots  in  Figures  7  and  8  where  it  is  clear  that 
the  M u  and  A/22  elements  are  more  similar  to  each  other  than  they  are  to  the  A/33  ele¬ 
ments  and  that  the  A/33  elements  have  additional  low  frequency  energy. 
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3.2.  Isotropic  Moment  Tensor 

The  primary  concern  of  this  paper  is  the  isotropic  part  of  the  moment  tensor, 
which  is  defined  here  in  terms  of  the  trace  of  the  moment  tensor  (Equation  2).  As 
mentioned  in  the  previous  section,  this  is  the  part  of  the  source  which  is  related  to  a 
change  of  volume  and  it  is  the  only  part  of  the  source  which  is  nonzero  for  an  ideal 
explosion.  The  isotropic  sources  which  were  estimated  for  the  two  events  Harzer  and 
Chancellor  are  shown  in  both  the  frequency  and  time  domain  in  Figure  9.  The  results 
are  quite  similar  for  the  two  events.  In  the  time  domain  these  isotropic  parts  are  fairly 
simple,  showing  slightly  acausal  one-sided  pulses  that  rise  to  a  peak  within  1  sec  and 
then  fall  to  a  constant  nonzero  level.  In  the  spectral  domain  they  show  high-frequency 
slopes  in  the  -2  to  -3  range,  spectral  comers  at  about  2  Hz,  and  additional  peaking  in 
the  specta  in  the  0.3-0.4  Hz  range.  This  peaking  in  the  spectra  below  1  Hz  tends  to 
obscure  the  comer  frequency  at  higher  frequency,  particularly  in  the  case  of  Harzer. 
The  interpretation  of  these  results  is  facilitated  by  comparing  them  with  the  predictions 
of  various  theoretical  and  empirical  models  of  explosive  sources. 

In  dealing  with  explosive  sources  it  is  convenient  to  consider  a  scalar  potential 
M'(t)  which  is  called  the  reduced  displacement  potential.  In  a  homogeneous  elastic 
medium  the  vector  displacement  u  at  any  time  t  and  position  x  is  simply 


u(f,x)  =  -  V  — - - 

x 


where  a  is  the  velocity  of  P  waves.  Because  of  this  simplicity,  much  of  the  previous 
work  on  explosive  sources  has  been  expressed  in  terms  of  the  reduced  displacement 
potential.  There  also  exists  the  following  simple  relationship  between  the  reduced  dis¬ 
placement  potential  and  the  isotropic  part  of  the  moment  tensor 

M,{t)  =  47tpa2vF(r)  (5) 

where  p  is  density.  Thus  the  estimates  of  the  isotropic  parts  of  the  moment  tensor 


TO 


SC* 5 


which  have  been  obtained  in  this  study  can  also  be  interpreted  as  estimates  of  the 
reduced  displacement  potential. 

The  reduced  displacement  potential  is  often  described  by  parameterized  models 
and  several  different  parameterizations  have  been  suggested.  Mueller  and  Murphy 
(1971)  considered  the  analytical  solution  for  a  pressure  pulse  acting  on  the  interior  of  a 
spherical  cavity  and  introduced  the  elastic  radius  of  the  source  and  the  pressure  acting 
at  that  radius  as  parameters.  The  characteristic  time  for  this  source  is  the  elastic  radius 
divided  by  the  velocity  of  compressional  waves.  Haskell  (1967)  introduced  an  empiri¬ 
cal  relationship  of  the  form 

H/(t)  =  TM[l  -/(f/x)]//(r)  (6) 

with 

fit)  =  e~‘  1  +  t  +  c2t2  +  c3t3  +  c4r4]  (7) 

and  where  x  is  a  characteristic  time  of  the  source,  H'00  is  the  long-time  response,  and  H 
is  a  unit  step  function.  Haskell  proposed  that  c2=  1/2,  c3=  1/6,  and  c4  be  an  adjustable 
parameter,  von  Seggem  and  Blandford  (1972)  suggested  an  expression  of  the  same 
form  but  with  c3=c4=0  and  c2  an  adjustable  parameter.  Helmberger  and  Hadley 
(1981)  also  adopted  the  same  form  and  proposed  that  c2=l/2,  c4=0,  and  c3  be  an  adju¬ 
stable  parameter.  Thus  we  have  four  different  functions  available  for  describing  the 
reduced  displacement  potential  and  each  functional  has  3  adjustable  parameters.  In  all 
cases  two  of  these  parameters  are  the  characteristic  time  x  and  the  long-time  level  4/00. 
For  the  Mueller-Murphy  model  the  third  parameter  is  the  ratio  of  peak  pressure  to  con¬ 
stant  pressure  at  the  elastic  radius,  while  for  the  other  three  models  is  is  one  of  the 
constants  in  equation  7. 

The  isotropic  parts  of  the  estimated  moment  tensors  shown  in  Figure  9  were  each 
interpreted  in  terms  of  the  four  different  models  of  the  reduced  displacement  potential 
that  are  described  in  the  previous  paragraph.  The  various  models  were  fit  to  the 
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observational  data  by  adjusting  the  parameters  of  the  model  until  what  appeared  to  the 
eye  as  a  reasonable  fit  was  achieved  in  the  time  domain.  In  essence,  the  characteristic 
time  x  was  determined  by  fitting  the  rise  time  of  the  first  pulse,  was  determined  by 
fitting  the  amplitude  of  the  first  pulse,  and  the  third  parameter  was  determined  by 
fitting  the  ratio  of  the  peak  amplitude  to  the  amplitude  in  the  interval  of  4-5  sec.  The 
results  of  this  fitting  process  are  shown  in  both  the  time  and  frequency  domain  in  Fig¬ 
ures  10-13  and  the  parameters  used  to  obtain  these  fits  are  given  in  Table  4. 

For  frequencies  less  that  1  Hz,  all  four  models  for  the  reduced  displacement 
potential  can  be  tit  to  the  isotropic  part  of  the  moment  tensor  so  as  to  satisfactorily 
explain  the  major  features  of  the  isotropic  moment  tensors  of  Harzer  and  Chancellor 
(Figures  10-13).  For  frequencies  above  1  Hz,  there  are  distinct  differences  between 
*he  abilities  of  the  models  to  fit  the  data.  The  Harzer  data  can  be  fit  better  than  the 
Chancellor  data  in  both  the  time  and  frequency  domains.  The  Mueller-Murphy  model 
does  slightly  better  than  the  other  models  in  fitting  the  Harzer  time  function  (Figure 
10),  particularly  for  the  back  side  of  the  first  pulse  between  1  and  2  sec.  In  the  fre¬ 
quency  domain  both  the  Mueller-Murphy  and  von  Seggem-Blandford  models,  which 
have  high-frequency  slopes  of  -2,  fit  the  Harzer  data  quite  well,  while  the  Haskell 
model  with  a  -4  slope  and  the  Helmberger-Hadley  model  with  a  -3  slope  underestimate 
the  high  frequency  energy  in  the  observational  data.  For  Chancellor  there  are  short- 
period  fluctuations  in  the  estimated  time  function  between  1  and  4  sec  which  vaguely 
suggest  the  possibility  of  a  multiple  event.  None  of  the  models  can  fit  these  fluctua¬ 
tions  and  all  of  the  models  underestimate  the  high-frequency  energy  in  the  Chancellor 
observed  spectrum,  particularly  in  the  1  to  5  Hz  range.  The  high-frequency  slope  of 
the  Chancellor  spectrum  is  between  -2  and  -3,  which  can  be  explained  by  the 
Mueller-Murphy,  von  Seggem-Blandford,  and  Helmberger-Hadley  models  but  not  by 
the  Haskell  model  Recalling  the  earlier  suggestion  that  both  a  low-frequency  function 
and  a  high-frequency  function  may  be  present  in  the  isotropic  part  of  the  moment 
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tensors,  it  is  apparent  in  Figures  10-13  that  the  process  of  fitting  the  models  to  obser¬ 
vational  data  in  the  time  domain  has  given  preference  to  the  low-frequency  function. 
This  is  particularly  evident  in  the  Chancellor  spectral  plots,  where  the  high-frequency 
function  with  a  comer  frequency  near  2  Hz  is  poorly  fit  by  the  models. 

Turning  now  to  the  model  parameters  which  were  estimated  during  the  fitting 
process  (Table  4),  we  see  that  the  values  of  these  parameters  can  depend  upon  the 
model.  The  estimates  for  vFao  and  x  show  variations  of  about  10%  and  30%,  respec¬ 
tively,  between  the  four  models.  The  pattern  is  consistent  for  x,  with  models  having 
steeper  high-frequency  slopes  requiring  shorter  characteristic  times.  This  model  depen¬ 
dent  variation  of  30%  in  x  could  lead  to  significant  differences  if  x  were  used  to  infer 
an  elastic  radius.  In  a  relative  sense  all  four  models  show  a  similar  pattern  in  that 
Harzer  has  a  which  is  about  25%  greater  than  for  Chancellor  and  a  x  which  is 
about  50%  greater. 

3.3.  Deviatoric  Moment  Tensor 

The  moment  tensor  of  an  ideal  explosion  would  have  only  an  isotropic  part  with 
all  deviatoric  parts  equal  to  zero.  This  is  not  the  case  for  either  of  the  events  Harzer 
or  Chancellor,  as  can  be  readily  observed  in  Figures  5,  7,  and  8  where  diagonal  ele¬ 
ments  are  not  all  identical  and  off-diagonal  elements  are  not  all  zero.  Even  more  basic 
evidence  concerning  the  existence  of  the  deviatoric  elements  can  be  inferred  from  the 
strong  early-arriving  energy  which  is  present  on  the  transverse  components  in  Figures 
2  and  3,  because  an  ideal  explosive  source  in  a  plane-layered  isotropic  medium  should 
generate  no  transverse  components.  In  fact,  noting  that  the  transverse  components  of 
ground  motion  are  comparable  in  amplitude  to  the  radial  and  vertical,  it  is  surprising 
that  the  off-diagonal  elements  of  the  estimated  moment  tensors  (Figures  7  and  8)  are 
so  much  smaller  than  the  diagonal  components.  What  appears  to  be  happening  is  that, 
either  because  they  are  incompatible  with  the  radial  and  vertical  components  or 


because  they  are  generated  by  a  mechanism  which  is  not  incorporated  in  the  calcula¬ 
tion  of  the  Green  functions,  the  transverse  components  of  ground  motion  have  a 
disproportionately  small  effect  upon  the  estimated  moment  tensors.  In  support  of  this 
interpretation,  note  that  in  Figure  6  the  observed  transverse  seismograms  are  con¬ 
sistently  underestimated  by  the  seismograms  calculated  with  the  estimated  moment  ten¬ 
sors. 

The  deviatoric  moment  tensor  (Equation  3)  is  a  symmetric  second-rank  tensor  and 
is  conveniently  interpreted  in  terms  of  its  3  eigenvalues  and  the  directions  of  the 
corresponding  eigenvectors.  Corresponding  to  the  usual  interpretation  of  principal  axes 
of  fault-plane  solutions,  these  eigenvalues  and  eigenvectors  have  been  identified  with 
the  tensional  principal  axis  T,  intermediate  principal  axis  I,  and  pressure  principal  axis 
P.  This  eigenvalue  decomposition  of  the  deviatoric  moment  tensor  as  a  function  of 
time  is  shown  in  Figure  14  for  both  Harzer  and  Chancellor. 

The  tensional  and  pressure  deviatoric  eigenvalues  in  Figure  14  are  only  slightly 
smaller  in  amplitude  than  the  isotropic  elements  in  Figure  9.  In  the  case  of  Harzer  the 
intermediate  eigenvalue  is  close  to  zero  throughout  the  first  5  sec  so  the  mechanism  for 
this  part  of  the  source  is  consistent  with  a  simple  shear  dislocation.  This  also  means 
that  the  tensional  and  pressure  eigenvalues  are  mirror  images  of  each  other  with 
respect  to  their  time  histories.  This  time  dependence  of  the  tensional  and  pressure 
eigenvalues  shows  a  general  similarity  to  the  isotropic  element  in  Figure  9.  Closer 
inspection  reveals  that  the  deviatoric  eigenvalues  are  delayed  in  time  by  about  0.5  sec 
with  respect  to  the  isotropic  element  and  have  longer  rise  times  of  about  1.5  sec  for 
the  first  pulse  as  compared  to  about  1.0  sec  on  the  isotropic  element.  In  the  case  of 
Chancellor  the  intermediate  eigenvalue  departs  from  zero  after  about  1  sec,  which  may 
indicate  a  source  mechanism  more  complicated  that  a  simple  shear  dislocation,  and  this 
causes  the  tensional  and  pressure  eigenvalues  to  differ  somewhat  in  their  time  his¬ 
tories.  However,  there  is  still  a  general  similarity  between  the  tensional  and  pressure 
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eigenvalues  and  the  isotropic  element  in  Figure  9.  As  in  the  case  of  Harzer,  the  devia- 
toric  eigenvalues  are  delayed  by  about  0.5  sec  with  respect  to  the  isotropic  element 
and  have  slightly  longer  rise  times  for  the  first  pulse. 

Referring  to  Figures  7  and  8,  one  can  see  that  the  size  and  time  histories  of  the 
deviatoric  eigenvalues  arises  primarily  from  the  differences  between  the  diagonal  ele¬ 
ments  of  the  moment  tensor  rather  than  the  amplitudes  of  the  off-diagonal  elements. 
Further  inspection  reveals  that  for  both  Harzer  and  Chancellor  the  major  cause  of  the 
deviatoric  parts  of  the  moment  tensors  is  the  fact  that  starting  about  0.5  sec  after  the 
origin  time  the  M33  elements  begin  to  differ  significantly  from  the  M  n  and  M 22  ele¬ 
ments.  This  explains  why  the  the  first  large  pulse  on  the  deviatoric  eigenvalues  is 
delayed  with  respect  to  the  first  large  pulse  on  the  isotropic  elements. 

While  Figure  14  describes  the  eigenvalues  of  the  deviatoric  parts  of  the  moment 
tensors.  Figure  15  describes  the  eigenvectors.  These  are  displayed  by  plotting  on  a 
stereographic  net  the  directions  of  the  eigenvectors  as  a  function  of  time.  As  these 
directions  can  be  quite  unstable  when  the  eigenvalues  are  small,  the  eigenvectors  are 
plotted  only  for  the  time  interval  0.5-3.0  sec  when  the  eigenvalues  have  large  values 
and  should  be  least  affected  by  noise.  A  simple  physical  interpretation  of  the  eigen¬ 
vector  directions  in  Figure  15  is  not  obvious.  Initially  the  P  axis  is  near  the  horizontal 
which  could  be  interpreted  as  strike-slip  motion  on  a  fault  plane  striking  northwest, 
and  this  roughly  corresponds  to  the  types  of  tectonic  stress  release  that  others  have 
found  for  Pahute  Mesa  explosions  (Wallace,  et  al.  1983,  1985;  Lay  et  at.,  1984).  As 
time  progresses  however,  the  P  axis  rotates  toward  the  vertical,  which  could  indicate  a 
transition  from  predominantly  strike-slip  to  predominantly  dip-slip  motion.  Given  the 
lack  of  an  obvious  simple  pattern  in  the  direction  of  the  eigenvectors,  the  fact  that  the 
deviatoric  elements  result  primarily  from  differences  between  the  larger  diagonal  ele¬ 
ments  of  the  moment  tensors,  and  the  fact  that  these  elements  have  a  spectral  content 
which  is  shifted  toward  the  low  frequencies  where  the  signal-to-noise  ratio  is  worse. 


considerable  uncertainty  must  be  attached  to  the  interpretation  of  the  deviatoric  parts  of 
these  moment  tensors. 


4.  INTERPRETATION 

The  moment  tensor  of  an  ideal  explosion  has  only  diagonal  elements  and  these 
three  elements  all  have  the  same  time  function.  The  moment  tensors  which  have  been 
estimated  for  Harzer  and  Chancellor  are  dominated  by  the  diagonal  elements,  but  they 
do  not  all  have  identical  time  functions.  In  particular,  the  M  33  elements  are 
significantly  different  from  the  Mn  and  M 22  elements.  Assuming  that  the  explosion 
began  as  an  isotropic  source  of  a  pressure  pulse  in  an  approximately  spherical  cavity, 
the  conclusion  follows  that  the  medium  has  responded  to  this  pressure  pulse  in  an 
asymmetric  manner,  with  the  response  in  the  vertical  direction  being  different  from 
that  in  the  horizontal  directions.  There  are  a  number  of  reasons  why  such  an  asym¬ 
metric  effect  might  be  expected  for  events  such  as  Harzer  and  Chancellor,  such  as  the 
vertical  effects  of  gravity,  the  strong  vertical  gradient  in  material  properties,  and  the 
fact  that  discontinuities  in  material  properties  are  predominantly  horizontal,  which 
includes  the  free  surface.  A  basic  question  is  whether  these  asymmetric  properties  of 
the  medium  are  sufficient  to  explain  the  asymmetric  effects  which  are  observed  in  the 
moment  tensor. 

An  effect  which  could  possibly  explain  the  asymmetry  of  the  Harzer  and  Chan¬ 
cellor  moment  tensors  is  the  phenomenon  known  as  spall.  For  buried  explosions  this 
effect  is  caused  primarily  by  the  direct  compressional  wave  being  reflected  from  the 
free  surface  with  a  change  in  sign,  which  can  result  in  tensional  stresses  at  particular 
depths,  nonlinear  failure,  part  of  the  material  separating  in  ballistic  flight,  and  a 
delayed  impact  when  this  material  rejoins  the  rest  of  the  earth.  Day  et  al.  (1983)  have 
modeled  spall  by  a  system  of  vertical  forces  in  which  momentum  is  conserved.  Stump 
(1985)  has  elaborated  upon  this  model  and  has  shown  that  for  high-frequency 
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waveforms  observed  near  the  source  the  contribution  from  the  spall  can  be  comparable 
in  amplitude  to  that  from  the  explosion.  The  spall  model  of  Day  et  al.  (1983)  consists 
entirely  of  vertical  forces  which  require  different  Green  functions  from  the  ones  used 
for  the  force  couples  of  the  second-order  moment  tensor  elements.  However,  Stump 
(1988)  has  shown  by  forward  modeling  that  the  Green  functions  for  the  vertical  force 
are  practically  indistinguishable  from  a  combination  of  the  Green  functions  of  the  diag¬ 
onal  elements  of  the  moment  tensor.  If  e(t)  is  the  time  function  of  the  explosion  and 
s(r)  is  the  integrated  time  function  of  the  forces  associated  with  the  spall,  then  Stump 
(1987)  suggests  that 

Mn  =  e(t )  +  s(t) 

M22  =  e(t)  +  s(t)  (8) 

Af33  =  £(r)  +  2s(r) 

If  spall  is  modelled  as  a  horizontal  tension  crack,  then  expressions  similar  to  equation 

8  are  obtained  except  that  the  factor  2  in  the  M33  term  becomes  ,  which  has  a 

A. 

value  of  2.9  at  the  depth  of  the  explosions  and  slightly  higher  values  at  shallower 
depths.  Both  of  these  results  suggest  that  the  presence  of  spall  is  consistent  with  the 
observation  that  the  M33  elements  of  the  moment  tensors  are  different  from  the  Mn 
and  Af  22  elements. 

Spall  is  commonly  observed  to  occur  with  shallow  events  such  as  Harzer  and 
Chancellor,  and  for  the  case  of  Harzer  there  is  considerable  quantitative  information 
available  concerning  the  details  of  the  process.  Patton  (1985)  has  studied  the  Harzer 
spall  process  using  data  from  in  situ  measurements,  pictures  of  the  ground  surface,  and 
computer  simulations.  He  finds  that  spall  for  this  event  was  a  fairly  symmetrical  pro¬ 
cess,  beginning  at  a  depth  of  290  m  about  0.3  sec  after  the  origin  lime  of  the  explosion 
and  extending  out  to  a  radial  distance  of  about  1.3  km.  The  ballistic  period  at  ground 
zero  was  between  0.9  and  1.7  sec,  with  this  range  in  interval  most  likely  due  to 
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multiple  spall  in  the  surficial  layers.  Patton  calculates  the  momentum  of  the  spall  mass 
as  1.51017  dyne-sec  and  finds  that  this  accounts  for  between  25%  and  60%  of  the  total 
isotropic  moment  of  the  Harzer  event,  where  the  range  depends  upon  the  ballistic 
period  assumed  for  the  spall  and  the  calculations  were  performed  for  the  excitation  of 
Lg  waves  in  the  frequency  range  of  0.2  to  0.5  Hz. 

Referring  to  the  moment  tensor  which  was  estimated  for  Harzer  (Figure  7),  the 
difference  between  the  M  33  element  and  the  other  diagonal  elements  is  generally  con¬ 
sistent  with  the  spall  process  which  has  been  outlined  by  Patton  (1985).  The  anomaly 
on  the  M 33  element  is  clearly  present  by  0.5  sec  and  lasts  until  about  2.5  sec,  which 
corresponds  roughly  to  the  time  period  in  which  the  spall  process  operated,  particularly 
if  one  allows  for  the  fact  that  the  spall  source  may  have  had  dimensions  of  2-3  km. 
The  sign  of  the  anomaly  changes  from  positive  to  negative  slightly  after  1  sec,  which 
may  correspond  to  the  beginning  of  spall  closure.  The  size  of  the  anomaly  is  roughly 
50%  of  the  isotropic  part  of  the  total  source,  which  is  also  consistent  with  Patton’s 
estimate  of  the  spall  contribution.  The  same  general  interpretation  can  also  be  given 
for  the  estimated  moment  tensor  for  Chancellor,  although  in  this  case  the  effects  of 
spall  appear  to  be  somewhat  smaller  than  for  Harzer  and  the  corroborating  information 
concerning  the  spall  process  is  not  available.  Thus,  on  the  basis  of  the  data  which  are 
available,  spall  must  receive  serious  consideration  as  a  possible  explanation  of  the 
asymmetry  in  the  diagonal  elements  of  the  moment  tensors. 

One  result  of  this  study  has  been  an  estimate  of  reduced  displacement  potentials 
from  the  analysis  or  near-source  waveforms  (Figure  9).  It  is  of  interest  to  compare 
these  results  with  other  estimates  of  reduced  displacement  potentials,  most  of  which 
have  been  based  on  free-field  measurements  at  buried  sites  within  a  few  hundred 
meters  of  the  sources  (Rodean,  1971,  1981;  Mueller  and  Murphy,  1971).  A  method  of 
making  such  comparisons  is  to  use  the  parameters  of  Table  4  which  have  been 
obtained  by  fitting  various  models  to  the  reduced  displacement  potentials.  This  is 
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particularly  convenient  for  the  Mueller- Murphy  model  because  a  fairly  complete  set  of 
scaling  relationships  are  available  for  this  model  (Mueller  and  Murphy,  1971;  Murphy, 
1977).  Using  the  scaling  relationship  of  Murphy  (1977,  Figure  15)  the  'F00  estimated 
for  Harzer  and  Chancellor  are  appropriate  for  mb  values  of  5.25  and  5.15,  respectively, 
which  are  about  0.25  less  than  the  published  values  (Table  1).  Stated  in  another  way, 
the  scaling  relationship  and  the  published  mb  values  predict  XF00  values  of  4-81010  cm3 
for  Harzer  and  3-61010  cm 3  for  Chancellor,  which  are  about  2  to  3  times  larger  than 
the  estimated  values  (Table  4).  Thus,  relative  to  scaling  relationships  based  on  the 
published  mb  values,  the  reduced  displacement  potentials  estimated  for  Harzer  and 
Chancellor  seem  to  be  too  small  by  a  factor  of  at  least  2.  There  are  a  couple  of  possi¬ 
ble  explanations  for  this  discrepancy. 

First  note  that  the  scaling  relationships  of  Murphy  (1977)  are  based  primarily 
upon  explosions  in  saturated  tuff-rhyolite  materials  and  the  explosions  Harzer  and 
Chancellor,  while  they  were  in  the  same  type  of  material,  were  located  above  the  water 
table  (Table  1).  Murphy  (1981)  and  Wallace  and  Barker  (1987)  have  both  found  that 
explosions  detonated  in  dry  tuff  and  alluvium  have  coupling  efficiencies  that  are  two  to 
three  times  less  than  explosions  detonated  in  saturated  tuff.  In  the  case  of  Murphy 
(1981)  this  difference  manifests  itself  in  mb  values  which  are  smaller  than  the 
saturated  tuff  predictions  for  a  given  yield  and  in  the  case  of  Wallace  and  Barker 
(1987)  it  manifests  itself  in  estimated  vF0O  values  which  are  smaller  than  the  saturated 
tuff  predictions  for  a  given  yield.  The  yields  of  Harzer  and  Chancellor  have  not  been 
announced  so  a  direct  comparison  with  these  published  results  is  not  possible,  but 
there  is  a  general  agreement  if  one  assumes  that  the  effect  of  dry  material  at  the  shot 
depth  is  to  decrease  'Poc  even  more  than  mb . 

Secondly,  it  is  important  to  recall  that  the  parameters  of  Table  4  were  determined 
by  fitting  the  isotropic  part  of  the  moment  tensor,  which  is  only  part  of  the  complete 
momen’  tensor  which  was  estimated.  The  mb  values  represent  the  effect  of  the 


complete  moment  tensor  which  is  considerably  more  than  just  the  isotropic  part  which 
is  represented  in  the  reduced  displacement  potential.  Additional  information  is  pro¬ 
vided  by  Patton  (1985)  who  used  Rayleigh  and  Lg  waves  recorded  at  regional  dis¬ 
tances  to  estimate  an  isotropic  moment  of  17- 1022  dyne-cm  for  the  Harzer  event, 
and  went  on  to  conclude  that  7-13-1022  dyne-cm  of  this  was  due  to  the  explosion  and 
the  rest  was  contributed  by  spall.  These  results  can  be  compared  to  an  average  value 
of  7.4- 1022  dyne-cm  from  Table  4  for  M{  of  Harzer,  which  is  about  a  factor  two  less 
than  what  Patton  (1985)  obtained.  Note  that  the  peak  value  of  the  spectrum  for  the 
Harzer  isotropic  moment  has  a  value  of  about  20-10 22  dyne-cm  at  0.3  Hz  (Figure  9) 
but  the  models  all  have  a  peaked  spectrum  with  a  lower  value  of  4/ao  (Figures  10-13). 

Another  parameter  of  Table  4  which  can  be  compared  with  scaling  relationships 
available  in  the  literature  is  the  characteristic  time  x.  Using  the  scaling  relationships  of 
Murphy  (1977,  1981)  the  mb  values  of  Harzer  and  Chancellor  can  be  converted  to 
predicted  values  of  0.14  sec  and  0.12  sec,  respectively,  for  the  characteristic  time  x.  In 
Table  4  we  find  estimated  values  based  on  the  Mueller-Murphy  model  of  0.55  sec  and 
0.35  sec,  respectively,  which  are  about  three  times  larger  than  the  predicted  values. 
Once  again  it  should  be  noted  that  the  scaling  relationships  being  used  here  were 
developed  for  explosions  detonated  in  saturated  tuff-rhyolite,  so  these  differences  may 
indicate  that  different  scaling  relationships  are  required  for  detonations  in  dry  materi¬ 
als.  Alternatively,  it  may  mean  that  the  depth  dependence  of  the  Mueller-Murphy  scal¬ 
ing  relationships  is  not  appropriate  for  explosions  above  the  water  table.  In  a  carefully 
controlled  study  of  small  chemical  explosions  Flynn  and  Stump  (1988)  found  that  as 
the  depth  of  the  explosion  decreased  the  amount  of  high-frequency  energy  decreased 
more  than  predicted  by  the  Mueller-Murphy  scaling  relationships,  and  they  attribute 
this  effect  to  the  change  from  full  to  partial  containment.  The  large  characteristic 
times  observed  in  this  study  are  qualitatively  consistent  with  the  results  of  Flynn  and 
Stump. 
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On  the  basis  of  the  comparisons  with  published  scaling  relationships  it  appears 
that  the  reduced  displacement  potentials  estimated  in  the  present  study  have  vF0O  values 
that  are  a  factor  of  about  2  too  small  and  x  values  which  are  about  a  factor  3  too  large. 
This  may  mean  that  the  method  of  moment  tensor  inversion  used  in  this  study  has  sys¬ 
tematically  underestimated  both  the  low-frequency  amplitude  and  the  high-frequency 
contributions  to  the  explosive  part  of  the  source,  or  it  may  mean  that  the  published 
scaling  relationships  can  not  be  applied  to  explosions  above  the  water  table  without 
modifications.  A  partial  answer  to  this  question  is  supplied  by  the  results  in  Figure  16. 

Figure  16  shows  the  observed  seismograms  and  several  different  versions  of 
predicted  seismograms  at  the  closest  available  site,  station  CIO  at  a  distance  of  1.8  km 
from  Chancellor.  The  closest  station  was  selected  with  the  idea  that  it  should  have  the 
minimum  propagation  effects  and  thus  provide  the  most  direct  estimate  of  the  source 
time  function.  The  seismograms  predicted  by  the  complete  moment  tensor  are  shown 
in  Figure  16  and  the  agreement  is  quite  good,  although  the  predicted  amplitudes  are 
slightly  less  than  observed  for  the  first  arrival.  Also  shown  in  Figure  16  are  the 
predicted  seismograms  which  are  obtained  when  only  the  diagonal  elements  of  the 
moment  tensor  are  included  and  when  only  the  isotropic  part  of  the  moment  tensor  is 
included.  By  comparing  these  results  with  those  for  the  complete  moment  tensor,  one 
can  see  the  progressive  effects  of  removing  the  off-diagonal  elements  of  the  moment 
tensor  and  then  removing  the  asymmetrical  part  of  the  diagonal  elements.  It  is  clear 
from  this  comparison  that  both  of  these  nonisotropic  effects  in  the  moment  tensor 
make  significant  contributions  to  the  predicted  seismograms,  including  the  first  arrival. 
Also  shown  in  Figure  16  are  the  seismograms  predicted  by  the  Mueller-Murphy  model 
of  an  explosion  using  the  parameters  of  Table  4.  For  this  model,  which  has  a  charac¬ 
teristic  time  of  0.35  sec,  the  predicted  amplitudes  are  smaller  than  the  observed  ampli¬ 
tudes  by  a  factor  of  2-3  for  the  first  arrival  and  the  first  downswing  is  not  well 
modeled.  Also  shown  are  the  predicted  seismograms  for  a  Mueller-Murphy  model 


25 


with  a  characteristic  time  of  0.18  sec,  which  improves  the  agreement  for  both  the 
amplitude  and  shape  of  the  first  arrival.  Consideration  of  observed  and  predicted 
seismograms  at  greater  distances  leads  to  the  same  general  result,  that  a  characteristic 
time  of  0.18  sec  predicts  the  first  arrival  on  the  observed  seismograms  better  than  the 
value  of  0.35  which  was  obtained  by  fitting  the  reduced  displacement  potential.  Simi¬ 
lar  results  were  obtained  for  the  other  models  of  the  reduced  displacement  potential 
which  are  listed  in  Table  4.  This  leads  to  the  conclusion  that,  if  a  simple  explosion 
were  assumed  for  the  source  and  only  the  first  cycle  on  the  radial  and  vertical  com¬ 
ponents  of  ground  motion  were  fit,  then  characteristic  times  would  be  obtained  which 
would  be  only  about  half  those  listed  in  Table  4,  which  bring  them  into  closer  agree¬ 
ment  with  the  published  scaling  relationships.  On  the  other  hand,  if  a  complete 
moment  tensor  is  used  as  the  source  model,  one  can  obtain  even  better  fits  of  the  first 
arrival,  do  a  much  better  job  of  fitting  the  later  arrivals  on  the  seismogram,  and  also  fit 
the  transverse  components. 

This  discrepancy  in  characteristic  times  is  related  to  the  discussion  of  section  3.2 
concerning  the  possibility  that  the  moment  tensors  contain  two  distinct  time  functions, 
a  low-frequency  function  which  is  strongest  on  the  A/33  element  and  which  controls 
the  modeling  of  the  reduced  displacement  potentials,  and  a  high-frequency  function 
with  a  comer  frequency  near  2  Hz  which  was  largely  ignored  in  the  fitting  process.  It 
now  appears  that  the  first  arrivals  on  the  seismograms  are  more  closely  related  to  the 
high-frequency  function  than  to  the  low-frequency  function.  This  could  be  interpreted 
to  mean  that  the  high-frequency  function  is  associated  with  the  explosion  while  the 
low-ffequency  function  is  associated  with  the  spall  process. 

Putting  all  of  the  these  results  together,  it  is  possible  to  describe  a  series  of  physi¬ 
cal  processes  which,  while  certainly  not  unique,  at  least  provide  one  plausible  interpre¬ 
tation  of  the  moment  tensors  which  were  estimated  for  Harzer  and  Chancellor.  The 
events  were  detonated  above  the  water  table  in  unsaturated  materials  and  spall  was  an 


important  part  of  the  process,  so  that  in  effect  these  were  only  partially  contained 
explosions.  The  spall  effect  contributed  an  asymmetry  to  the  source,  beginning  at  the 
latest  a  few  tenths  of  a  second  after  detonation,  so  that  the  effective  forces  in  the  verti¬ 
cal  direction  had  different  amplitudes  and  different  time  functions  than  those  in  the 
horizontal  directions.  The  moment  tensor  estimates  were  constructed  on  the  basis  of 
the  radiated  seismic  waveforms,  and  because  none  of  the  stations  were  at  distances  less 
than  three  source  depths,  all  the  processes  in  the  immediate  source  area,  including  both 
the  explosion  and  spall,  were  interpreted  as  being  caused  by  a  single  effective  source 
process.  Thus  the  moment  tensors  contained  the  combined  effects  of  these  various 
processes,  including  the  explosion  which  may  have  had  characteristic  times  of  as  small 
as  0.1  sec  and  the  spall  which  may  have  had  characteristic  times  of  as  long  as  1  sec  or 
more.  Extracting  the  isotropic  part  of  the  moment  tensor  helps  isolate  the  explosive 
part  of  the  source,  but  it  is  not  completely  successful  because  the  spall  process 
involves  a  change  in  volume  and  thus  also  appears  in  the  isotropic  part.  It  appears 
that  the  reduced  displacement  potentials  which  have  been  estimated  in  this  study  con¬ 
tain  contributions  from  spall  in  addition  to  the  explosive  source,  which  results  in  a 
possible  overestimation  of  the  characteristic  time  of  the  explosion.  The  inelastic  spall 
process  absorbs  energy  which  might  have  otherwise  gone  into  seisnuc  waves,  resulting 
in  a  decreased  seismic  coupling  coefficient  and  a  reduced  displacement  potential  which 
is  smaller  than  what  is  usually  obtained  for  fully  contained  explosions  in  saturated 
materials. 
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S.  DISCUSSION 

The  problem  of  deciphering  the  physics  of  explosive  sources  from  the  study  of 
radiated  elastic  waves  has  been  separated  into  two  different  inverse  problems  in  this 
study.  The  first  inverse  problem  is  to  remove  the  propagation  effects  and  obtain  an 
estimate  of  the  effective  force  system  operating  in  the  immediate  vicinity  of  the 
source,  which  has  been  characterized  by  the  second-order  moment  tensor.  Unique 
solutions  are  possible  for  this  inverse  problem,  and  in  the  case  of  the  explosions 
Harzer  and  Chancellor  stable  results  were  obtained  in  the  frequency  range  of  0.2  to  5.0 
Hz.  This  was  possible  because  of  the  large  number  of  well-distributed  recording  sites 
and  the  availability  of  a  good  estimate  for  the  average  velocity  structure  in  the  region 
of  study.  The  second  inverse  problem  consists  of  interpreting  the  moment  tensor  in 
terms  of  the  physical  processes  acting  in  the  source  region.  Solutions  to  this  problem 
are  inherently  nonunique,  so  additional  information  and  assumptions  must  generally  be 
included  in  the  interpretation  process. 

A  basic  result  for  the  explosions  Harzer  and  Chancellor  is  that  the  moment  tensor 
is  dominated  by  its  diagonal  elements  but  that  these  elements  contain  an  important 
asymmetry,  namely  that  the  M33  element  differs  significantly  from  the  Mu  and  M 22 
elements.  Considerable  evidence  points  to  the  possibility  that  this  asymmetry  is 
related  to  spall  processes  taking  place  in  the  depth  interval  between  the  explosion  and 
the  free  surface.  An  interpretation  in  terms  of  a  source  which  contains  the  effects  of 
both  a  symmetric  explosion  and  an  asymmetrical  spall  is  capable  of  quantitatively 
explaining  many  of  the  features  of  the  moment  tensor.  For  distances  in  the  range  of  1 
to  11  km  and  frequencies  in  the  range  of  0.2  to  5.0  Hz,  these  two  different  source 
effects  have  comparable  amplitudes  and  it  is  not  easy  to  separate  their  relative  contri¬ 
butions  to  the  moment  tensor.  Thus,  if  the  isotropic  part  of  the  moment  tensor  is  inter¬ 
preted  as  being  due  solely  to  a  simple  explosion,  the  results  could  be  misleading. 
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The  tentative  conclusion  that  both  the  explosion  and  spall  have  made  significant 
contributions  to  the  radiated  elastic  waves  which  were  analyzed  in  this  study  calls  into 
question  one  of  the  basic  assumptions  of  the  method  which  was  used,  namely  that  the 
source  is  sufficiently  localized  in  space  so  that  it  can  be  represented  by  a  second-order 
moment  tensor.  It  would  appear  that  an  inversion  should  be  performed  which  takes 
into  account  the  fact  that  two  sources  may  be  acting  to  two  different  locations,  one  at 
the  detonation  point  of  the  explosion  and  the  other  at  a  shallower  depth  where  spall  is 
suspected.  While  this  is  possible.  Stump  (1988)  has  already  performed  numerical 
experiments  for  this  type  of  inverse  problem  and  the  results  are  not  promising.  He 
finds  that  when  the  closest  stations  are  several  times  the  separation  of  the  two  sources, 
there  is  a  strong  trade-off  between  the  two  sources  and  the  inversion  in  highly 
nonunique.  Having  observations  from  points  closer  to  the  sources  might  improve  the 
situation  but  one  then  encounters  another  problem,  that  of  interpreting  ground  motions 
in  the  nonlinear  zone  of  the  source  where  the  basic  wave  equations  which  have  been 
assumed  throughout  this  study  are  no  longer  valid. 

This  paper  has  concentrated  on  the  interpretation  of  the  most  prominent  parts  of 
the  moment  tensor,  the  diagonal  elements  The  deviatoric  elements  were  also  analyzed 
but  a  definitive  interpretation  could  not  be  given.  The  events  Harzer  and  Chancellor 
were  both  detonated  in  parts  of  Pahute  Mesa  where  there  have  been  numerous  previous 
explosions  and  neither  seemed  to  show  evidence  of  strong  tectonic  release.  It  would 
be  helpful  if  the  same  type  of  analysis  used  in  this  study  could  be  applied  to  an  explo¬ 
sion  known  to  have  significant  tectonic  release.  It  could  be  that,  given  the  experience 
gained  in  studying  the  moment  tensor  from  such  an  event,  one  could  then  return  and 
give  a  more  meaningful  interpretation  of  the  deviatoric  elements  of  the  Harzer  and 
Chancellor  explosions. 

Finally,  it  should  be  pointed  out  that,  while  the  transverse  components  of  ground 
motion  were  included  in  the  inversion,  helped  constrain  the  estimates  of  the  moment 


tensors,  and  were  partially  modeled  by  the  predicted  seismograms,  they  have  not  been 
satisfactorily  explained.  In  fact,  the  interpretation  in  terms  of  a  symmetric  explosion 
and  a  vertical  spall  process  predicts  no  transverse  motion,  a  result  obviously  at  vari¬ 
ance  with  the  observations  of  transverse  motions  comparable  in  amplitude  to  the  radial 
and  vertical  motions  at  all  distance  ranges.  Apparently,  our  understanding  of  the  gen¬ 
eration  of  elastic  waves  by  explosions  is  still  incomplete. 
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Table  1.  Shot  Parameters 


Harzer 

Chancellor 

Date 

6  June  1981 

1  Sept  1983 

Origin  Time 

18h  00m  00.1s 

14h  00m  00.1s 

Latitude 

37.303N 

37.273N 

Longitude 

1 16.326W 

1 16.355W 

Surface  elevation 

2097m 

2040m 

Depth  to  shot 

637m 

625m 

Shot  Medium 

tuff 

tuff 

Depth  to  standing  water 

668m 

647m 

Magnitude  ML  (BRK) 

5.4 

5.3 

Magnitude  mb  (USGS) 

5.5 

5.4 

Table  3.  Network  Data  for  Chancellor 

$ 

Station  Distance 

Azimuth 

Elevation 

First  Peak  Acceleration 

Maximum  Acceleration 

(km) 

(deg  E  of  N) 

(m) 

(cm/sec  ) 

(cm/ see  } 

•  Vj 

R  T  V 

R  T  V 

5 

Table  4.  Parameters  for  Fitting  Reduced  Displacement  Potentials 


Harzer  Chancellor 


von  Seggem-Blandford 

Mj  „  (1020  dyne-cm) 

740 

570 

(101(W) 

2.7 

2.1 

X  (sec) 

0.50 

0.35 

c2 

-2.5 

-1.5 

Hclmbcrger- Hadley 


M/ao  (1020  dyne-cm) 

(1010cm3) 

X  (sec) 
c2 
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Figure  2. 
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Figure  4. 


Figure  5. 


Figure  6. 


Figure  7. 


Figure  8. 


Figure  9. 
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Map  showing  the  locations  of  the  explosions  and  recording  sites 
in  the  Pahute  Mesa  region  of  the  Nevada  Test  Site.  The  two 
explosions  Harzer  and  Chancellor  are  denoted  by  open  circles 
and  the  sites  of  temporary  recording  stations  are  denoted  by 
closed  circles. 

Velocities  of  ground  motion  for  the  event  Harzer  which  were 
obtained  by  integrating  the  observed  accelerations.  The  plots 
have  been  scaled  by  multiplying  by  the  epicentral  distance 
(Table  2).  The  radial  direction  R  is  positive  away  from  the 
source,  the  transverse  direction  T  is  clockwise  about  the  source 
as  viewed  from  above,  and  the  vertical  direction  V  is  positive 
upward.  The  seismograms  have  been  shifted  to  align  the  first 
arrivals  which  corresponds  to  a  reducing  velocity  of  about  4 
km/sec. 

Similar  to  Figure  2  for  the  event  Chancellor. 

Estimates  of  average  material  properties  for  the  upper  4  km  of 
the  Silent  Canyon  Caldera. 

Amplitude  density  spectra  of  the  moment  rate  tensors  estimated 
for  the  events  Harzer  and  Chancellor.  The  solid  lines  are  for  the 
Mn  elements  of  the  moment  tensors,  the  dashed  lines  for  the 
M33  dements,  and  the  dotted  lines  for  the  average  deviatoric  ele¬ 
ments.  In  the  lower  parts  of  the  plots  are  shown  the  estimated 
standard  errors  which  were  obtained  by  assuming  random  white 
noise  in  the  input  data. 

Comparison  between  the  observed  and  predicted  ground  veloci¬ 
ties  for  two  stations  of  the  event  Chancellor.  For  each  com¬ 
ponent  the  predicted  result  is  shown  directly  below  the  observed 
data.  The  amplitude  scales  are  different  for  each  component  but 
are  the  same  for  observed  and  predicted  traces.  The  panel  on  the 
left  is  for  station  C3  at  an  epicentral  distance  of  2.1  km  and  that 
on  the  right  is  for  station  C7  at  a  distance  of  9.8  km. 


Time  domain  estimates  of  the  6  elements  of  the  second-order 
moment  tensor  for  the  event  Harzer. 


Time  domain  estimates  of  the  6  elements  of  the  second-order 
moment  tensor  for  the  event  Chancellor. 


The  spectral  amplitude  density  of  the  isotropic  parts  of  the 
moment-rate  tensor  and  the  time  domain  plots  of  the  isotropic 
parts  of  the  moment  tensor  for  the  explosions  Harzer  and 


w 


Figure  10. 

Figure  11. 

Figure  12. 

Figure  13. 

Figure  14. 
Figure  15. 

Figure  16. 


Chancellor. 

Amplitude  and  spectral  amplitude  density  (solid  lines)  of  the  iso¬ 
tropic  part  of  the  estimated  moment  tensors  and  the  Mueller- 
Murphy  model  of  the  reduced  displacement  potential  (dashed 
lines)which  were  fit  to  the  moment  tensor  results  for  the  events 
Harzer  and  Chancellor.  The  parameters  used  to  achieve  these 
fits  are  listed  in  Table  4. 

Amplitude  and  spectral  amplitude  density  (solid  lines)  of  the  iso¬ 
tropic  part  of  the  estimated  moment  tensors  and  the  von 
Seggem-Blandford  model  of  the  reduced  displacement  potential 
(dashed  lines)which  were  fit  to  the  moment  tensor  results  for  the 
events  Harzer  and  Chancellor.  The  parameters  used  to  achieve 
these  fits  are  listed  in  Table  4. 

Amplitude  and  spectral  amplitude  density  (solid  lines)  of  the  iso¬ 
tropic  part  of  the  estimated  moment  tensors  and  the  Helmberger- 
Hadley  model  of  the  reduced  displacement  potential  (dashed 
lines)which  were  fit  to  the  moment  tensor  results  for  the  events 
Harzer  and  Chancellor.  The  parameters  used  to  achieve  these 
fits  are  listed  in  Table  4. 

Amplitude  and  spectral  amplitude  density  (solid  lines)  of  the  iso¬ 
tropic  part  of  the  estimated  moment  tensors  and  the  Haskell 
model  of  the  reduced  displacement  potential  (dashed  lines)which 
were  fit  to  the  moment  tensor  results  for  the  events  Harzer  and 
Chancellor.  The  parameters  used  to  achieve  these  fits  are  listed 
in  Table  4. 

Eigenvalues  of  the  deviatoric  elements  of  the  estimated  moment 
tensors  for  the  events  Harzer  and  Chancellor. 

Lower-hemisphere  equal-area  stereographic  plot  of  the  directions 
of  the  eigenvectors  of  the  estimated  deviatoric  moment  tensor  for 
the  events  Harzer  and  Chancellor.  The  directions  are  shown  for 
the  time  interval  0.5-3  sec.  Also  plotted  as  triangles  are  the 
positions  of  the  recording  sites  as  they  are  projected  onto  the 
focal  sphere  along  the  path  of  the  direct  P  wave. 

Observed  and  predicted  seismograms  at  station  CIO  at  a  distance 
of  1.8  km  from  the  explosion  Chancellor. 
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